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ABSTRACT 

We outline how redshift-space distortions (RSD) can be measured from the angular correla- 
tion function w{0), of galaxies selected from photometric surveys. The natural degeneracy 
between RSD and galaxy bias can be minimized by comparing results from bins with top-hat 
galaxy selection in redshift, and bins based on the radial position of galaxy pair centres. This 
comparison can also be used to test the accuracy of the photometric redshifts. The presence of 
RSD will be clearly detectable with the next generation of photometric redshift surveys. We 
show that the Dark Energy Survey (DES) will be able to measure f(z)a$(z) to a ler accuracy 
of (17 x b)%, using galaxies drawn from a single narrow redshift slice centered at z = 1. 
Here b is the linear bias, and / is the logarithmic rate of change of the linear growth rate with 
respect to the scale factor. Extending to measurements of w(8) for a series of bins of width 
0.02(1 + z) over 0.5 < z < 1.4 will measure 7 to a la accuracy of 25%, given the model 
/ = £l m (z) 7 , and assuming a linear bias model that evolves such that b = 0.5 + z (and fixing 
other cosmological parameters). The accuracy of our analytic predictions is confirmed using 
mock catalogs drawn from simulations conducted by the MICE collaboration. 

Key words: Cosmology - Observational: methods, large-scale structure 



1 INTRODUCTION 

The statistical analysis of redshift-space distortions (RSD) provides 
constraints on the build-up of cosmological structure (see H amilton] 
|1998| for a review). RSD are apparent anisotropic patterns of galax- 
ies formed when redshifts are converted into galaxy distances as- 
suming the Hubble flow: coherent comoving peculiar velocities are 
translated into apparent distributions of galaxies. On small distance 
scales (~ 1 h -1 Mpc), we see the well-known "finger of god" ef- 
fect (see, e.g., |Jackson| 1972} , due to orbital velocities of galaxies 
that are members of bound structures. On larger distance scales, 
galaxies move with structure growth, leading to an apparent excess 
in the clustering strength along the line-of-sight. 

Following linear theory, the large-scale effect of RSD on 2- 
point statistics can be modelled as a simple enhancement of the 
galaxy overdensity fie ld jKaiser 1987) . This leads to an anisotropic 
correlation function (Hamilton 19921 that is a function of both 
radial and transverse separations. The excess line-of-sight clus- 
tering strength depends on the amplitude of the matter veloc- 
ity field that, following linear theory, is usually parametrised by 
/(2)o"8 !maS s(z), where / is the rate of change of the linear growth 
rate / = d log G/d log a, G is the linear growth rate, a is the scale 
factor of the Universe, and as,-ma B a{z) is the root mean square am- 
plitude of fluctuations in the matter overdensity field in spheres of 
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radius 8 ft -1 Mpc. We will now drop the "mass" label as we will 
only refer to ag (z) as normalising the matter fluctuation field from 
here-on. Without the RSD, the amplitude of the real-space galaxy 
overdensity field can be parametrised by b(z)as(z), where 6(2) is 
the bias of the galaxies. For ACDM cosmologies, to a good ap- 
proximation, f(z) = f2 m (2) 7 iLinder 2005 1 with 7 = 0.557, and 



measurements of £ s (r) can be used to test the standard cosmologi- 
cal model — either by translating them into measurements of f(z) 
using an external constraint on the galaxy bias or performing a dual 
fit for b(z)ag(z) and f(z)crg(z) (while accounting for the implied 
change in cr&{z) when f{z) is changed). 

RSD measurements have traditionally been confined to spec- 
troscopic redshift surveys, where accurate redshifts have been ob- 
tained for a sample of galaxies: measurements from the 2dF Galaxy 
Redshift Survey (2dFG RS;|Colless et al.|2003] > and the Sloan Dig- 
ital Sky S urvey (SDSS;|Ybrk et al.|2000[> aregi ven by |Peacock et 
|aTJ(200T);|Hawkins et al.H2003);|Percival et al.|l|2004");|Pope et al. 



l|2004);|Zehavi et al.| ( |2005| l; |Okumura et al.H2008> ; |Cabre & Gaz 
tanaga ( 2009 ) for example, and RSD have recently been detected at 
higher redshift ( |Guzzo et al.|2009)|Blake et al.|2010| >. 

Photometric surveys are usually overlooked when it comes 
to RSD measurements because the photometric redshift errors 
are thought to wash out any RSD signal. Without question, 
these errors will significantly degrade the imprint of RSD. How- 
ever, it is timely to investigate the extent to which this is 
true given both the high-quality of upcoming data from the 
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Dark Energy Survey (DES www.darkenergysurvey.org), 
the Panoramic Survey Telescope and Rapid Response System 
(PanStarrs pan-starrs.ifa.hawaii.edu), and the Large 
Synoptic Survey Telescope (LSST www . Isst . org), and recent 
theoretical analyses that have investigated the effects of RSD on 
projected clustering measurements. 

If galaxy redshifts are used to select samples, then we 
see redshift-space effects in angular/projected clustering measure- 
ments, even though the angular positions of galaxies are unchanged 
( |Fisher et al.fl993| [Nock et al.|2010| hereafter N10). Corrections 
for these effects were recently applied to angular power-spectra 
analyses of the SPSS by |Blake et al.| < |2007| >, |Padmanabhan et al.| 
l |2007) , and |Thomas et al.| ( |2010| >; thelatter of whom have used the 
corrections to place constraints on /. 

Working in configuration space, these corrections have been 
shown to agree with w(9) and its covariance calculated with mock 
photometric redshift catalogs (N10, Cro cce et al.||2010b| here- 
after CIO). This effect arises because RSD impart coherent fluc- 
tuations in galaxy density along any redshift-dependent boundary 
imposed on the galaxies. Correct treatment of the redshift-space 
distortions is extremely important, as studies suggest that the most 
precise (non-RSD) cosmological constraints are obtained from us- 
ing narrow redshift bins ( Sanch ez~et al.|2010| >, which also give the 
strongest RSD signal (N10, CIO). 

N10 have shown that constructing galaxy samples using pair- 
centre binning removes the effects of RSD. Such a binning scheme 
does not consider the individual redshift of a galaxy, but rather the 
average redshift of (or distance to) the pair of galaxies. Thus, a pair- 
centre binning scheme with 0.45 < z < 0.55 will include all of the 
galaxy pairs with average redshifts between 0.45 and 0.55. In stan- 
dard 'top-hat' binning, all of the galaxies with redshifts between 
0.45 and 0.55 would be included. The symmetry imparted by the 
pair centre binning is sufficient to remove the effect of RSD. 

RSD and distortions caused by using the wrong geometri- 
cal model when analysing data in three dimensions (the Alcock- 
Palczynski (AP) effect Alcock & Paczynski 1979) give rise to sim- 
ilar distortions in the measured clustering (Ballinger et al. 1996 
|Simpson"ife Peacock 2009 ), although the strength of the degeneracy 
depends on the cosmological assumptions used in the data analysis 
( |Samushia et al.pOll) . Thus, in general, analyses of RSD using 
3D data either have to be performed simultaneously with geometri- 
cal measurements, or utilise (and rely on) more precise geometrical 
constraints. The measured angular correlation function, which only 
depends on redshifts and angles, does not depend on the cosmo- 
logical model to be fitted. Fitting models to the measured angular 
clustering is therefore considerably simpler. 

In this work we investigate the precision to which one can 
measure cosmological structure growth using galaxy samples with 
photometric redshifts. In Section|2] we present the theoretical mod- 
elling we employ to estimate angular correlation functions of galax- 
ies, w(9) and their error, for different galaxy binning schemes. In 
Section [3] we test our analytic error predictions against mock cat- 
alogs produced by the MICE simulation and show the importance 
of including non-linear effects when measuring f(z)as(z). In Sec- 
tion]^ we determine the precision to which f(z)as(z) can be mea- 
sured for a survey similar to DES using redshift bins of different 
sizes and mean redshifts and galaxies with various mean photomet- 
ric redshift errors and linear bias. In Section[5] we test the robust- 
ness of the measurements against changes in the cosmology and 
systematics in the photometric redshifts and we show how compar- 
ing measurements using standard binning schemes and those using 
pair-centre binning measurements can test the accuracy of the pho- 



tometric redshifts. We conclude in Section [6] Where appropriate 
we assume a flat ACDM model with Q m = 0.25, h = 0.7, and 
n b = 0.045. 



2 MODELING ANGULAR CORRELATION FUNCTIONS 

For a sample of galaxy pairs, the angular-correlation function 
w(8), can be modeled by integrating the redshift-space correlation- 
function £ s (r), over the joint distribution of galaxy pairs /(zi, Z2), 
( |Phillippsetal.|1978^ : 



(6) = / dz\ I dz 2 f(zi, z 2 )£ s [r(6, z\, z 2 ) 



(1) 



where the galaxy separation r is a function of the angular separa- 
tion of the galaxies 6 and their redshifts z\ and zi. For a binning 
scheme based on galaxy position (such as top-hat binning), we can 
separate the selection function for the two galaxies. We can then re- 
place f(zi, 22) by n{z\)n{z 2 ), where n{z) is the normalized true 
redshift distribution, in the above expression. We note that, in prac- 
tice, n{z) is an estimate usually obtained from the photometric red- 
shifts and their errors. We investigate potential systematics due to 
this estimation in Section [5TT1 

For the pair-centre binning, f(z\,z 2 ) can be replaced by 
n pc (zi)rid(\zi — 22I), where n pc (z) is the normalized distribu- 
tion of pair centres and nd(Az) is the normalized distribution of 
pair separations (in redshift; note that the tophat w(0) is recovered 
when integrating over n pc (zi)nd(\zi — z 2 \) when zi and z 2 are 
restricted to lie within the tophat bounds). We incorporate the pho- 
tometric redshifts into f(z\ , z 2 ), working with a redshift-space cor- 
relation function, rather than photometric redshift-space one, so the 
functions n(z) and f(z\,z 2 ) depend on the photometric redshift 
error. For simplicity these errors are assumed to be Gaussian with 
standard deviation a z . In practice, this means that we Monte Carlo 
sample the true redshift distribution twice, and then Monte Carlo 
sample the respective Gaussians to find two positions in photomet- 
ric redshift space. If the pair-centre of these positions lies within 
the photometric redshift bin, the true redshift pair-centre and true 
redshift separation is added to its respective distribution (and the 
respective distributions are then normalized). See N10 for further 
details. 

In this work, we only consider the plane-parallel approxima- 



tion to the geometry, and assume that £ s (r) is given by 1 Hamilton 
[T992l > 

?Qi,r)=Mr)PoiM) + b(r)&(M) + U(r)P*<Ji), (2) 
where 

Mr) = (& 2 + §&/ + i/ 2 )£M, 



6 to = (f&/+|/ a )fc(r)-£'(r)], 



(3) 
(4) 



&(r) = ^/ 2 K(r) + ^'(r)-^"(r)], (5) 
Pi are the standard Legendre polynomials, and 



£' = 3r" 3 / i{r'){r') 2 dr' 



(6) 



i" = 5r- J / ar')(r'fdr', (7) 
Jo 

b is the large-scale bias of the galaxy population being considered 
and £(r) is the isotropic 3-dimensional, real-space correlation func- 
tion. The factor /i is the cosine of the angle between the separation 
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along the line of sight and the transverse separation. We calculate 
£(r) by Fourier transform from a power spectrum calculated as de- 
scribed in Eisenstein & Hu ( 1998) and include the non-linear mod- 
elling given in equation 10 of CIO. 

The model we employ makes some strong assumptions (linear 
RSD, linear bias, weak non-linear gravity). We note that its appli- 
cation to projected/angular correlation functions is well tested with 
simulation (N 10, CIO) and has been shown to be a good fit to mea- 



sured w(8) (see, e.g., Ross et al. 2007 2008 2010b for discussions 



of the range of scales over which measurements of angular clus- 
tering are well-fit by a linear bias model and, e.g., Samushia et al. 
1201 II for a discussion of the limits to which the linear RSD model 
is valid). 

Each moment of £(r) entering Eq.[2]is normalised to the stan- 
dard correlation function, and hence its amplitude is fixed by <rf (2) . 
Any measurements of b(z) or f(z) calculated using this model are 
therefore perfectly degenerate with the value of erg (2). It is there- 
fore most convenient for our purposes to parametrize this model 



with two, separate, factors: b(z)as(z) and f(z)as(z) (e.g. Song & 
|Percival|200"9| >. 

Eq. [T] shows that binning schemes act as "moments" of the 
anisotropic redshift-space correlation function: each provides an in- 
tegral of this function weighted by a different distribution of angles 
to the line-of-sight. The binning schemes therefore work in a sim- 
ilar way to the Legendre polynomials in the standard definition of 
the anisotropic correlation function. Comparing results from differ- 
ent schemes such as considering measurements from both the top- 
hat and pair-centre schemes, would allow the simultaneous mea- 
surement of both b(z)ag(z) and f(z)ag(z) (when accounting for 
the covariance between these measurements). 



2.1 Errors and Covariance 

Recent work has considered the problem of analytically calculating 
covariances for angular/projected clustering measurements l |Cohn| 
|2009| C10). In order to set constraints on how well future surveys 
will measure the angular clustering, we follow the approach of C 10, 
which is based on the angular power spectrum for redshift slice i, 
which we denote Ci ,i. 

The error on the full-sky angular power spectrum, Ce, is given 
by (see, e.g. , |Dodelson|200"3] > 



(d) ~ {c t ){c t ) 

2 {C t + 1/nf 



(2^+1) 



(8) 
(9) 



where n is number density of objects per steradian. The predicted 
variance of the cross-power spectrum, C\' 3 , between samples i and 
j is given by (see, e.g., eqs. 18 and 19 of |White et al.) 2009 1 



(2«+l). . (in) 

x [( C ^) 2 + (Cj + l/nO(^ + l/^)] . 



Thus, given that 



t>0 



V 4tt ) 



Pi(cos6)C; 



(11) 



and that the effect of only observing a fraction of the sky, fsky, is 
well approximated by simply dividing a 2 (Ce) by f B k y (see, e.g., 
C10), the covariance (that we denote Cov) between measurements 
of the angular cross-correlation function at 8 and 6' between red- 
shift slices i and j (that is the auto-correlation in the case i = j) 



can be modelled as 

cov^, e , 9 , = {wi^wijie')) - {w id (9))(wi,j(e')) (12) 

= E*>o S^Pi(cos6)P e (co S e')a 2 (C^) + (13) 

where n v is the number of galaxy pairs (that can be determined 
from the number density and angular size of the bin) and 5g t gi is 
the Kronecker delta. 

Given P(k), the Ce^j spectra can be calculated via (see, e.g., 
|Fisher et al.|1993||Pa7jnianabhan et al.|2007) 



c e j = ~ J fc 2 dfcP(fc)* f , l (fc)*,, J (fc), 
where \f^ ( i is given by 

9 e ,i(k) = ^li + J dzni(z)D(z)jt(kx(z)), 

and je is the spherical Bessel function of order £ and 



*I,i(*) = //6 / dzn t (z)D(z) 



(2l+3)(2t-l) 
(l+l)(l+2) 



(21-1) (2£+l) ^- 2 " (2<+lj(2l+3) ^+2( fc X( Z )) 



(14) 



(15) 



(16) 



This set of equations provides the full covariance matrix in both 
scale and redshift, as required in order to make measurements of 
f(z)cr&(z) at multiple redshifts. 



2.2 Pair-Centre Covariance 

For the pair-centre binning, determining the covariance is not as 
straightforward. We can produce an approximate solution by cal- 
culating the variances of all of the thin slices in photometric red- 
shift that would contribute to the pair-centre measurement and all 
of the co- variances between these slices. Given a photometric red- 
shift survey, one can split it into very narrow consecutive shells of 
width Ss in photometric redshift located at 



Ss x i (i — 0, 00). 



(17) 



If one wishes to estimate the pair-centre angular correlation func- 
tion, Wpc, in a bin centered at Sb of width Ast (in photometric 
redshift space), the condition for a pair of galaxies lying in bins i 
and j to enter the selection is 



Vij = \ s ij — s b \/As b < 1 



(18) 



with Sij = Ss x (i + j)/2 the mean photo-z of the galaxy pair. 
Thus the condition can be implemented by a 0(y) function giving 
1 if y ^ 1 and otherwise. 



w pc (8) 



The pair-centre correlation can be estimated as, 



E*E,-;>i n ( 2 »M 2 j)©(jft: 



(19) 



where Wij(6) is the angular cross-correlation function between 
bins i, j and n(z,) is the value of the overall redshift distribu- 
tion n(z) (i.e., the redshift distribution of sources over the entire 
survey being considered) at redshift Zi with C\ 3 calculated as de- 
scribed by Eqs. 14-16, with rii(z) in Eq. 15 being the true redshift 
distribution of the thin slice in photometric redshift space, and the 



covariance of Wij(8), Covij (#,#'), is therefore given by Eq. 13 
Thus, one can calculate all of the GaVij(9, 6') contributing to the 
In order to determine Covd C (6I 8'), one must ac 



19 



sum in Eq. 

count for the covariance between Wh,i(8) and Wj : k{8), which we 
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o h 



(20) 



denote CoVh,i,j,k- We can estimate Co"Vh,i,j,k by calculating the 
"overlap", Oh,i,j,k, of the four thin photometric redshift shells, 

J dzn h (z)ni(z)nj(z)nk(z) 
J dzn 2 h (z)n 2 (z) 

and assuming 

Cov h ,i,j,k = Oh,i,j,kCov h ,i(n = oo) (21) 

where we set n = oo since the covariance between the separate Ce 
measurements should not have a shot-noise contribution. We thus 
estimate 



ii". 



CoVp 

where we use . . as shorthand for . . and 



Wt = n{z h )n{z i )n{z j )n{z k )Q(y hi )Q{y jk ). 



(23) 



This approximation is validated by comparisons with mock cata- 
logs in section[3] 

We estimate the covariance between the top hat 
and pair-centre binning by assuming the correlation, 
Cov pc> TH / ^CovpcCovTH, between the measurements is 
equal to the ratio of the effective volume (given by Eq. 5 of 
|Tegmark]|1997} of the galaxies used in each binning scheme. We 
find this correlation varies between 1/7 and 1/2, depending on the 
specific bin widths and median redshift. We find that these values 
are consistent with the covariance we find in the mocks used in 
Section |3~Tj though the mocks do exhibit some scale dependence. 
Allowing the covariance to change by as much 25% does not affect 
any of the conclusions of this work, but this issue is worthy of 
further study in the future. 



2.3 Estimating Errors for DES 

As a specific example, we now consider the sensitivity of the an- 
gular correlation functions expected to be measured for DES. As 
in N10, we take the overall selection of DES galaxies to have the 
form 



riDEs(z) oc ( ( fg) 2ex P( _ (^) 



(24) 



and assume that DES will observe 300 million galaxies over l/8th 
of the sk)Q These galaxies are assumed to be unbiased on average, 
that is b — 1, We will consider the dependence of our conclusions 
on bias in Section l4~T1 

Because the DES galaxy density is high, the predicted error 
for clustering measurements at a given equivalent physical scale 
(i.e, ~ x( 2 )#> where xi z ) ls the comoving distance to redshift z), 
for bins that are of the same order as a z , will be driven almost en- 
tirely by the median redshift of the sample: the shot noise term in 
Eq. [T3]has a negligible impact on the final error. Thus, for DES, 
Cavg gi oc Cf , and we would expect the same percentage error 
regardless of the width of the redshift bin for both top-hat and pair- 
centre binning. Note that this argument will break down around 
the BAO scale where the angular correlation crosses from positive 
to negative and consequently the proportionality arguments do not 



1 We thank the DES LSS working group for providing these estimations. 



(22) 
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Figure 1. Model angular correlations functions (multiplied by 9 2 for clar- 
ity), calculated directly using Eq. l[TJ from the 3D correlation function. 
These are for distributions of galaxies within the photometric redshift bin 
0.45 < z < 0.55 and assume a z = 0.03(1 + z). The black line is for a 
top-hat bin when redshift-distortions (RSD) are included, its and theoretical 
uncertainties are also displayed (calculated via Eq. |13) , while the red line 
was calculated without RSD. The solid blue line is for a pair-centre binning, 
with its theoretical uncertainties (calculated via Eq. |22) also displayed. The 
dashed-blue line displays the top-hat, no RSD prediction divided by 5.62. 



hold. Also, moving to smaller bins will not necessarily provide sig- 
nificant additional signal because results will be correlated through 
the photometric redshift errors. 

We assume Gaussian photometric redshift errors, a z , of 
0.03(1 + z). Although the average o z of the full sample of DES 
galaxies at these redshifts is predicted to be closer to 0.05(1 + z) 
than 0.03(1 + z) (see |Banerji et al.|2008| Fig. 1 of C10), a slight 
reduction in the number of galaxies used can have a large effect on 
the mean value of a z . As discussed above, the density of galaxies 
is high so results are cosmic variance and photometric redshift lim- 
ited: we can therefore reduce the sample density without a signifi- 
cant increase in noise. But the ability to measure RSD does depend 
strongly on the width of the redshift distribution: Consequently, 
galaxy samples should be constructed to have minimal average o z 
(by selecting, for example, primarily early-type galaxies). The ease 
of redshift determination for these galaxies can lead to significant 
gains: SDSS galaxies with 0.2 < z < 0.3 and M r < -21.2 have 
an average a z of 0.041, but if one removes all of the galaxies with 
a z greater than 0.053(1 + 2), this is reduced to 0.03 and 82% of the 
galaxies are left in the sample. (This sample was studied by |Ross et] 
|al.|2010a| and the estimated photometric redshift errors were con- 
sistent with the measured clustering.) 

Fig. [T] compares model correlation functions (multiplied by 
9 2 for clarity), calculated by directly projecting the 3D correlation 
function. These were calculated for galaxies within a photometric 
redshift bin 0.45 < 2 < 0.55 and a z = 0.03(1 + 2). We plot mod- 
els for a top-hat bin in real-space (red) and redshift-space (black), 
and for pair-centre binning (blue). We also display the theoretical 
uncertainties, as predicted via Eq.[T3]for the top-hat bin and Eq.|22| 
for the pair-centre bin. The signal-to-noise is considerably worse 
for the pair-centre binning. 
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0(degrees) 

Figure 2. Model angular correlation functions for galaxies drawn from the 
DES in top-hat bins with photometric redshifts 1.0 < z p hot < 1-13 
(black) and 0.45 < z p hot < 0.55 (red). The solid curves include RSD, 
dashed lines do not. Error bars, which are only plotted for models with 
RSD, were calculated using Eq. j!3( . 



The dashed line in Fig.[T]represents the case where we divide 
the real-space top-hat prediction by 5.62. This illustrates that the 
shape of the pair-centre measurement is nearly identical to that of 
the top-hat without redshift space distortions. We find this to be 
true in general when the width of the pair-centre bin matches the 
width of the top-hat bin. While we find we can apply slightly more 
narrow redshift bounds to the pair-centre bin and obtain models that 
are slightly better matched, for simplicity, we will always apply the 
same redshift bounds when comparing model pair-centre and top- 
hat w(0). 

Fig.[2]presents the predicted angular correlation functions for 
DES, with and without RSD, for top-hat photometric redshift bins 
centred at z — 0.5 and z = 1.065. For clarity, we only present 
errors for the models with RSD. Data are plotted for two redshift 
bins, and we clearly see that the increase in volume surveyed at 
high redshifts leads to more precise clustering measurements. 



3 TESTING WITH MOCKS 

This section outlines tests of our analytic prediction for the model 
w(6) and its covariance matrices, performed using mock catalogs 
drawn from the N-body simulation MICE7680 produced by the 
MICE collaboration. This simulation assumed a flat ACDM cos- 
mology with parameters fi m = 0.25, fib = 0.044, ag jmass (0) = 
0.8, n 3 — 0.95 and h = 0.7 and tracked the gravitational evo- 
lution of 2048 3 dark-matter particles within a como ving volume 



Fosalba et al. 



with a box size given by Lbox = 7680 h 1 Mpc (see 
[2008l and |Crocce et al.|2010a| for further details). We use mock cat- 
alogues corresponding to dark matter particles distributed in slices 
of varying width centered at the mean redshift of z=0.5 (i.e. ex- 
tracted from the co-moving output at this z) and subtending 1/8 of 
the sky. The very large volume of MICE7680 allow us to obtain at 



Theory 
Mocks 



1 2 3 4 5 

0(degrees) 

Figure 3. The variance (a 2 ) of the angular correlation function in the pair- 
centre bin 0.475 < z pc < 0.525 for objects with a DES-like distribution, 
determined from 125 mock realizations (black triangles) and estimated the- 
oretically (using Eq.|22j . 



least 125 independent mocks for every bin- width and photometric 
error under consideratiorjf] We refer the reader to C10 for a detailed 
presentation of the mock build-up. 



3.1 Pair-centre Covariance 

We use the mocks to test our estimation of the pair-centre covari- 
ance matrix. We note that the theoretical modelling of wth and its 
covariance matrix has been extensively tested against the mocks in 
C10. Using a range in redshift 0.2 < z < 0.8, the volume was 
split, in photometric-redshift-space, into shells 15/i _1 Mpc thick. 
We then calculated w pc (8) for each realization by measuring all 
of the correlation- and cross-correlation functions of shells where 

0. 475 < (zi+Zj)/2 < 0.525, where m is the mean redshift of slice 

1, and then summing them according to the overall redshift distribu- 
tion. These measurements allow us to calculate the covariance ma- 
trix for w pc (9) for the pair-centre bin 0.475 < z pc < 0.525. Fig.[3] 
presents the diagonal elements of the covariance matrix, multiplied 
by 10 s , determined from the mocks (red squares) and from Eq.|22| 
using shells 0.01 thick in photometric-redshift-space (black trian- 
gles). The theory appears to slightly underestimate the variance on 
large scales and slightly over-estimate the variance on small scales. 
However, the agreement is close enough that we should be able to 
use Eq.|22]for the purposes of this study. We find that the agreement 
between off-diagonal terms is quite similar, in terms of the absolute 
differences between the theoretical and mock estimates. 



3.2 Constraining f(z) ag(z) 

While the model wth and its covariance matrix has already been 
shown to agree with the mocks (see C10), we wish to confirm 



Mock catalogs are publicly available at http://www.ice.cat/mice 
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Figure 4. Histograms displaying the distribution of /(z)<r 8 (z) values ob- 
tained from 125 MICE mock catalogs with z = 0.5 and Az = 0.15(l + z) 
(solid black line) and Az = 0.05(1 + z) (dashed red line) with no 
redshift error, for a mock catalog with Az = 0.03(1 + z) and red- 
shift errors 0.02(1 + z) (dotted blue line), and for a mock catalog with 
Az = 0.05(1 + z) and redshift errors 0.02(1 + z) (long-dashed green 
line). 



that we can recover an unbiased value of f(z)as(z) and its un- 
certainty. We have tested four mock samples against the analytic 
predictions: two that assume no photometric redshift errors, with 
widths of Az = 0.05(1 + z) and Az = 0.15(1 + z) and two 
that have redshift errors 0.02(1 + z) and Az = 0.03(1 + z) and 
Az = 0.05(1 + z) (all have RSD effects applied). We only test 
the predictions at medium redshift (z ~ 0.5), where we can select 
125 non-overlapping samples from the MICE simulations. While 
we are testing the results at z = 0.5, we expect the DES to obtain 
much more precise constraints at higher redshifts. However, given 
that the distribution of matter becomes more non-Gaussian as the 
Universe evolves, we can reasonably expect that the modelling that 
works at z ~ 0.5 should work at higher redshift as well. Thus, by 
confirming that our modelling is correct at z ~ 0.5, we expect our 
modelling to be correct in general. 

The mocks provide measurements of w(8) with 0.2° ^ 
8 ^ 8° in angular bins of width 0.2°. Jack-knife errors, 
<jjk(8), (see, e.g., Scranton et al. 2002 Myers et al. 2007 



|Ross et al.|[2007> are calculated for each of these 125 mocks 
using 63 independent regions of equal area. In order to esti- 
mate the covariance matrix of each realization, we determine 
the covariance matrix of the ensemble of mocks, Cov e , and 
use the covariance matrix, Covjfc given by Covjk(8,8') — 

Cov e (6,6')crj k (0)cTj k ((0')/y/Cov e (0,6)Cov e (8',0'). We find 
that, in the case of the most extreme fluctuations from the mean, 
the size of the jack-knife errors of individual angular bins are cor- 
related with the size of the fluctuation. Modulating the covariance 
matrix based on the jack-knife errors thus reduces the weight of the 
more extreme fluctuations and we find this yields results matching 
our analytic predictions. 

To find the best-fit values of f{z)as{z), we fit for data in 



Figure 5. Histograms displaying the distribution of /(z)<r 8 (z) values ob- 
tained from MICE mock catalogs with z = 0.5, Az = 0.05(1 + z), and 
o z = 0.02(1 + z), when non-linear effects are included in the modelling 
(solid black) and are not included in the modelling (dashed red). 



the range 1° < 8 < 5°, marginalizing over 6(2)0-8(2). In or- 
der to produce the models we compare to these mocks, we use 
the same power-spectrum used by the MICE collaboration (avail- 
able at the MICE website). Fig. [4] presents histograms of these 
best-fit f(z)a$(z) values obtained for four types of mock; Az = 
0.15(1 + 2) in solid black, Az = 0.05(1 + 2) in dashed red, 
Az = 0.03(1 + 2) & a z = 0.02(1 + 2) in dotted blue, and 
Az = 0.05(1 + 2) & a z = 0.02(1 + 2) in long dashed green. 
For the mean and standard deviation of the two distributions with- 
out redshift errors, we find 0.47±0.09 for Az = 0.05(1 + 2) 
and 0.48,±0.15 for Az = 0.15(1 + 2). For a z = 0.02(1 + 2), 
we find 0.42±0.18 for Az = 0.03(1 + 2) and 0.43±0.2 for 
A z = 0.05(1 + 2). 

We make analytic predictions, based on the modelling de- 
scribed in Section [2] using methods described fully in Sec- 
tion HI These predictions are /(0.5)o- 8 (0.5) = 0.44 ± 0.09, 



/(0.5)a 8 (0.5) 



0.44 



+0.13 



/(0.5)a 8 (0.5) 



0.44+°;if, and 



/(0.5)<7 8 (0.5) = 0.441q 2o, respectively, for the four redshift bins 
described above. These agree quite well with the results from the 
mocks and therefore suggest that our analytic formalism can be 
used to effectively predict the precision to which f(z)ag(z) can be 
measured with a DES-like survey. We further note that we have ap- 
plied a linear RSD model (given by Eqs. 2-7) to obtain all of these 
results, and this does not appear to bias the results. 



3.3 Importance of Non-linear Gravitational effects 

The MICE simulations obviously include non-linear effects from 
the gravitational clustering in the distribution of matter, and we 
need to model these effects in order to match the recovered angu- 
lar correlation function. If we do not add non-linear effects to the 
spatial correlation function (i.e., £(r, 2) = D 2 (z)£n n (r, 0)) when 
fitting to the data, we obtain a significantly different distribution of 
best-fit f(z)ag(z) values. We note that in both cases we are still 
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Figure 6. The mean angular correlation function from the MICE mocks 
for a redshift bin centred on z = 0.5, Az = 0.05(1 + z), and a z = 
0.02(1 + z) is shown by the black error-bars. The solid black curve displays 
the model w(8) including non-linear effects and the LCDM prediction of 
/(0.5)<rg(0.5) = 0.44. The dotted red curve displays the model when non- 
linear effects are not included and /(0.5)<rg(0.5) = 0.44 and the dashed 
blue curve displays the model when non-linear effects are not included and 
/ (0.5) erg (0.5) = 0.37 (that is the best-fit when non-linear effects are not 
included). 



using the linear model described by Eqs. 2-7 to apply the effects 
of RSD. The histogram of this distribution for the redshift bin with 
Az = 0.05(1 + z) and a z — 0.02(1 + z) is displayed displayed 
in red in Fig. [5] (with the original histogram including non-linear 
effects displayed in black). The distribution is significantly skewed 
towards low f(z)a%(z) values. The mean of the distribution de- 
creases to 0.37 (from 0.43) — a 14 per cent change. 

Fig.|6]shows model w(8) x 8 2 (curves) and the meanw(9) x 8 2 
from the MICE mocks (black error-bars) for a redshift bin cen- 
tred at z = 0.5, Az(l + z) = 0.05, and a z = 0.02. The 
black curve displays the model including non-linear effects and 
/(0.5)<rg(0.5) = 0.44 (that is the standard ACDM prediction at 
z = 0.5). The red curve displays the model without non-linear ef- 
fects and /(0.5)<78(0.5) = 0.44, and the blue-curve displays the 
same model for /(0.5)<rg(0.5) = 0.37 (that is the best-fit when 
the non-linear effects are not included). Essentially, the non-linear 
effects soften the BAO peak, and when they are not included in the 
models, the data is best-fit by a model with a smaller f(z)as(z) in 
order to compensate. We note that the error-bars displayed are not 
the standard deviations, ctm, of the mocks but the errors in their 
mean (i.e., -t=ctm and there are N = 125 mock realizations). 

V N 

Thus, the disagreement at large scales between the model (black 
curve) and the mocks is at the ~ la level. Given that the f(z)a$(z) 
values determined from the mocks agree with the analytic predic- 
tions, this disagreement at large scales is too slight to significantly 
bias our results. 

We find that the best-fit value of /(0.5)cr 8 (0.5) has shifted by 
14 per cent when we do not include non-linear effects. The main 
non-linear contribution is due to the damping of the BAO due to 



Figure 7. The solid line displays the model f(z)ag(z) for our default 
ACDM model. 1<t errors (black) were calculated for the expected mea- 
surements made via successive top-hat photometric redshift bins for DES 
galaxies between 0.475 < z < 1.42 and the blue error-bars display Fisher 
matrix predictions for similar redshift bins. The red points with lcr errors 
are the measurements made with the WiggleZ survey 4Blakeetal.|2010) . 

large-scale velocity flows. This effect thus grows as the Universe 
evolves (see, e.g., Appendix A of C10). This implies that we expect 
the linear theory to yield a less biased result at higher redshift. This 
is indeed the case. We can use our theoretical predictions to find 
the linear model w(0) that is the best-fit to the non-linear model 
w(8). At z — 0.7, we find decreasing the value of f (0.7) as (0.7) 
by 8 per cent best matches the non-linear model, and that at z — 
1.0, decreasing the value of /(1.0)<rg(1.0) by 6 per cent yields 
the best-fit to the non-linear model. While the magnitude of the 
bias does decrease with redshift, it is quite clear that accurate non- 
linear modelling of w(8) is required to obtain accurate values of 
f(z)a 8 (z). 



4 MEASURING REDSHIFT-SPACE DISTORTIONS IN 
THE DES 

Fig. [2] shows that RSD will significantly affect the angular clus- 
tering measured from DES for top-hat bins. We now quantify the 
constraints that could be possible for the DES. The blue error-bars 
in Fig.|7]represent the uncertainty that one obtains from Fisher ma- 
trix predictions for the full-anisotropic 3D power spectrum for the 
DES. These were calculated using the method described in |White| 
et al. (2009), with a radial damping term for the power spectrum 
of the form e _ ' fc<Tzf1 ' , where we assume the photometric redshift 
error a z is the dominant cause of radial damping (and /i remains 
the cosine of the angle to the line of sight). These predictions are 
for successive bins of width 0.03(l+z), the minimum redshift of the 
first bin being 0.5, and a z = 0.03(1 + z). 

We now consider how closely we can match the Fisher ma- 
trix errors using angular clustering measurements. We shall use the 
angular correlation function for pair centre binning, w pc , which is 



© 2010 RAS, MNRAS OQO,[7|fT3| 



8 A. J. Ross et al. 



independent of RSD, to measure b{z)as(z). The angular correla- 
tion function for top-hat binning depends on both f(z)as(z) and 
b(z)as(z), so the combination of the two measures both. Thus, de- 
spite the relative weakness of the signal-to-noise of w pc (see Fig. 
[TJ, the fact that w pc is independent of f(z)as [z] allows us to break 
the f(z)cr&(z) - b(z)a$(z) degeneracy inherent in most clustering 
measures. (Though, it should be noted that the Fisher matrix pre- 
diction does not include the information about pairs outside of the 
photometric redshift slice.) 

We assume that all binned correlation function measurements 
are drawn from a multi-variate Gaussian distribution, with covari- 
ance matrices calculated as described in Section|2T] Thus, we cal- 
culate the standard \ 2 statistic, constructing the full covariance ma- 
trix from those of wth (Eq. [T3| and w pc (Eq. 22) and estimating 
Covth, P c as described at the end of section 2.2| We are there- 
fore able to find the range of f(z)as(z) values that take up 68 
per cent of the likelihood space around the Likelihood Maximum, 
while marginalising over the uncertainty to which the bias can be 
measured. 

We constrain the range of angles used in these constraints to 
be such that 20k" 1 Mpc< r eq < lOO/i" 1 Mpc, where r eq = 
2x(z)ta,n(9/2) (with 8 in radians) and x(z) is the co-moving dis- 
tance to the median redshift, z, of sample of galaxies we consider 
(this has been shown by, e.g., CIO to be an appropriately small 
scale to which the modelling is still accurate). We use 15 contigu- 
ous top-hat photometric redshift bins, and 15 matched pair-centre 
bins between 0.475 < z < 1.42, each of width 0.0333(1 + z), so 
the first bin has 0.475 < z < 0.525 and the last 1.34 < z < 1.42. 
(The fact that Fig. 1 of CIO suggests the mean photometric redshift 
uncertainties are expected to ~ constant over the range 0.45 < z < 
1.4 justifies our chosen redshift range.) We conservatively assume 
that we can select 7 million galaxies in each redshift bin, with good 
photometric redshifts (a z ~ 0.03(1 + z)) out of a total of 300 
million galaxies: the approximate DES galaxy density distribution 
shows that we should expect more than 10 million galaxies in each 
bin in total. 

The expected constraints from these angular clustering mea- 
surements are shown in Fig. [7] compared with the Fisher matrix 
predictions. At high redshift, the predictions are nearly identical, 
but at low redshift the Fisher matrix predicts better constraints than 
using the angular clustering method described above. This is con- 
sidered further in the next section, where we determine that the op- 
timal bin at z ~ 1 is ~ 0.03(1 + z), but at z ~ 0.5, a narrower bin 
recovers better constraints. We have also plotted recent measure- 
ments made using the WiggleZ survey (Blake et al. 2010) in Fig.|7] 
WiggleZ already does better than DES at its redshifts, but beyond 
the redshift limits of WiggleZ (z ~ 0.8), DES will provide compet- 
itive measurements with better than ~20% error for slices of width 
Az = 0.03(1 + z). We note that the f(z)a$(z) constraints be- 
tween separate redshift bins are highly correlated with each other. 
We take these correlations into account when determining how well 
the DES constrains f(z)as(z) over its full redshift range in section 

E2 



4.1 Optimizing the Method 

The size of the effect of RSD on angular clustering measurements 
is strongly dependent on the characteristics of the sample, and in 
particular the true radial size of the redshift slice is one of the most 
important factors. This is given by a convolution between the pho- 
tometric redshift bounds of the bin and the photometric redshift 
distributions. The mean bias of the galaxies is also important, as 
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Figure 8. The black points show the expected error on /(z)crg(z) versus 
the width of the photometric redshift bin, for unbiased tracers with average 
photometric redshift error cr z = 0.03(1 + z), selected from redshift bins 
centred on z = 1.065. The other lines display the same information, for 
samples changed as labeled. 



the relative amplitude (compared to the real-space clustering) of 
the redshift distortion effect is inversely proportional to the bias. 

Fig. [8] displays the expected error on f(z)ag(z), for various 
possible data samples. The black line is for unbiased tracers cen- 
tred at z = 1.065 and with a z — 0.03(1 + z), which gives a mini- 
mum in the expected error for the bin Az(l + z) — 0.033(1 + z). 
This represents a balance between the number of galaxies in each 
sample against the size of the RSD signal. The blue line shows the 
effect of a sample selected with larger photometric redshift errors, 
now 0.05(1 + z). The expected error on f(z)a$(z) is more than 
twice as high, and the minimum is at a larger bin width. The true 
bin width is determined by a convolution between the photometric 
redshift bounds of the bin and the photometric redshift distributions 
so decreasing the photometric redshift limits has less of an effect 
when (7 z is increased. 

Fig. [5] also shows that the effect of increasing b has a detri- 
mental effect on the RSD measurements: a 50 per cent increase in 
the average galaxy bias results in a ~ 50 per cent increase in the 
predicted RSD error. (The error on /3 = f/b remains constant, im- 
plying the error f(z)as(z) is proportional to the bias.) This is an 
important consideration for the selection of galaxy samples, partic- 
ularly at high redshifts where the samples have larger biases (see, 
e.g., |Zheng et al.|2007) , as do red galaxies that currently have the 
lowest photometric redshift errors. 

At lower redshifts, Fig. [8] shows that RSD measurements are 
strongly degraded by the significantly smaller volume available. At 
these redshifts, the error decreases significantly with bin size, down 
to the most narrow bin we test, Az(l + z) = 0.02(1 + z), though 
it does appear to asymptote to ~ 0.16. This is due to the fact that 
the number densities are so high at this redshift that decreasing the 
total number of galaxies has a minimal effect on the error on w{9). 
It is clear from Fig. [8] that the choosing the optimal sample with 
which to measure RSD depends on a complicated balance of many 
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interdependent factors: this can only reliably be achieved after all 
of the data is in hand. 



4.2 Cosmological Constraints from full Redshift Range 

We now consider cosmological constraints from the combination of 
measurements at different redshifts. We account for the covariance 
in f(z)a$(z) between redshift bins by assuming that their Pear- 
son coefficient is equal to the Pearson coefficient that defines the 
covariance of w(6) between different redshift bins. We determine 
the covariance of w(0) between different redshift bins by assuming 
(see, e.g., [Thomas et al. 2010) 

Cov(Ci,Ci) = j^ Y) (Ci' j f (25) 

and then substituting Cov(C| , C{ ) for a 2 (C^ ) in Eq.fo] yielding 
Cov(u>(#i)\ u>(#2)"'). We find that the Pearson coefficient given 
by Cov(w(6-l)\ w(6 2 y )/ Cov Mjei : g 2 ,e 2 is nearly in- 

dependent of 61,62- Thus, we set 

CoV/,j,j / J GoVf,i,iCoVf,j,j = Cov(lo\ W^)/y/ C0Vi,iC0Vjj(26) 

where Covfjj is the covariance in f(z)crg(z) between redshift 
bins i and j. We find that if we make Gaussian realizations of 
our simulated data matching the covariance matrix for correlation 
functions expected for two adjacent redshift bins, the covariance 
between measurements of f(z)as(z) recovered is within 10% of 
what we estimate using Eq. |26| (and is not systematically offset). 
Thus, we use Eq. [26] to construct our full covariance matrix for 
f(z)as(z), given its relative ease in computation and that it should 
be sufficiently accurate for our tests. 
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Figure 9. 68% confidence limits for a constant offset in f(z)crg(z) from 
our fiducial model, when we combine the constraints produced by a series 
of Wi(9) measured in successive redshift bins with median redshifts be- 
tween 0.5 and 1.4, versus the width of the photometric redshift bins. The 
black symbols show the prediction for unbiased tracers with average pho- 
tometric redshift error cr z = 0.03(1 + z). The red points display the same 
information for <r z = 0.05(1 + z) and the blue points assume that the bias 
of the selectable galaxies evolves, such that b(z) = 0.5 + z and the average 
photometric redshift error of these galaxies is cr z = 0.03(1 + z). 



4.2.1 Constant offset for f (z)ag(z) 

If we consider deviations around our fiducial model with a redshift- 
independent offset in f(z)as(z), then the 68% confidence interval 
for Af (z)os(z) is given by 



A(/(*)<7 B (*)) = 1/ 



(27) 



We determine A (f(z)ag(z)) for unbiased tracers occupying a 
variety of bin widths (all with the median redshift of the first 
bin at 0.5 and the median redshift of the final bin less than 1.4; 
the number of galaxies in each bin is assumed to be 7xl0 6 x 
(Az(l + z)/0.05)) and plot the results in Fig.|9] We find that, for 
cr z = 0.03(1 + z), A (f(z)os(z)) asymptotes towards a value of 
0.025. For a z = 0.05(1 + z) the decrease in the error as a function 
of bin width is much stronger, but at the minimum bin width we 
probe (0.02(1 + z)) we find A (f(z)o 8 (z)) = 0.06 — more than 
twice as large as for o z = 0.03(1 + z). 

We also investigate a scenario with more realistic evolution in 
the bias of the galaxies selected, assuming b(z) = 0.5 + 2. This al- 
lows the bias to be unity at z = 0.5 and gives a similar value to the 
bias of galaxies in the DEEP2 field (see, e.g., Zheng et al.|2007| > at 
z = 1. For this bias model, we find A (f(z)as{z)) asymptotes to- 
wards 0.04. Overall, this suggests that DES alone should be able to 
rule out models that cause the growth history to differ by more than 
10 per cent from the ACDM prediction. It further suggests that, 
though the covariance between bins grows large (the reduced co- 
variance of off-diagonal terms is as large as 0.9), extra information 
can be gained by making the bins as narrow as possible. 



4.2.2 Deviations in 7 

One can also determine how well DES will be able to constrain 
models of the type f(z) — fi m (z) 7 (see, e.g., Linder |2005) . We 
note that this also implies a change in the behaviour of 0$ (2), since 
as (z) = o 8 (0)G(z) and f(z) = dlog(G)/dlog(a). We fix the 
value of (7 8 (1088) that yields cr 8 (0) = 0.8 in the ACDM cosmol- 
ogy and determine the expected behaviour of f(z)as(z). We then 
test this against the f(z)os(z) covariance matrix we obtain for a 
series of Wi(9) measurements to predict the ability of DES to con- 
strain 7. The results we find are presented in Fig.[l0] We find there 
is a steady decrease in the expected error on 7 as the bin width 
deceases, and we find that DES alone should be able to place la- 
constraints on the value of 7 to less than 0.25, even with realistic 
evolution in the bias of the observed galaxies. 



5 ROBUSTNESS TO UNDERLYING COSMOLOGY AND 
PHOTOMETRIC REDSHIFT SYSTEMATICS 

The RSD constraints presented in the previous section assumed per- 
fect knowledge of cosmological geometry and photometric redshift 
errors. Both changing the assumed cosmology and introducing sys- 
tematic errors can cause scale dependent effects that could mimic 
the effects of redshift-space distortions when measuring w(0) using 
a top-hat bin. 

Measurements of the angular correlation function only de- 
pend on angles and redshifts, and are therefore independent of 
cosmological model: a model does not have to be specified in or- 
der to make these measurements. The cosmological dependence 
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Figure 10. As Fig.|5] but now showing the expected constraints on the value 
of 7, given the model f(z) = O m (z) 7 . 



only arises through the model to be fitted to the data. We test the 
systematic effect of assuming an incorrect background model on 
f(z)ag(z) measurements by using the 1 < z v hot < 1.13 redshift 
bin. In effect, we consider fitting to a fixed ACDM model with free 
structure growth, and determine the degree to which the inferred 
structure growth changes when the assumed ACDM parameters are 
incorrect. 

We calculate the value of f(z)a%{z) one measures if one as- 
sumes the geometrical cosmology given by Q m = 0.25, ft = 
fib/fim = 0.18, h — 0.7, but the Universe has a true cosmology 
that is different. For (fb)true ~ 0.2, the best-fit value of f(z)as(z) 
increases to 0.52±0.07 (from the true value of 0.42±0.07) while 
it decreases to f(z)a 8 (z) = 0.31± o ; o | for {f b )tr M = 0.16. 
Changing Q m produces an even more significant effect, as the 
f(z)ag(z) decreases to 0.231° io when (^m)i™ e = 0.3. While 
these changes to the cosmology are within the 95% confidence 
limits determined by WMAP (Komatsu et al. 2009), by the time 
that DES finishes, the Planck Experiment will have set tighter 
constraints on these parameters, making this effect less important. 
However, given current constraints, a joint fit to all of the parame- 
ters is clearly required. 

There are other effects that can cause scale dependent changes 
to the observed w(0). Any scale dependent bias will clearly be 
a problem, and it has recently been shown that large halos dis- 
play scale dependent biases (see, e.g., Smith et al. 2008, Manera| 



|& Gaztanaga||20"09"l |Des]acques et al.|2010) . However, as shown 
in Fig. 8, the best RSD constraints will be obtained from objects 
that reside in the least biased halos, for which we expect the least 
scale-dependent bias. Any non-zero }nl will also cause a scale 
dependent bias, though for |/jvl| ~ 100, this is a much smaller 
effect than RSD at scales smaller than the BAO scale. Additionally, 
we may worry about observational systematics, such as Galactic 
extinction or stellar contamination. Ross et al. ( 2011) has shown 
that a thorough accounting of and correction for such systematics 
and their effects is possible, but depending on the particular survey, 
this could translate to additional systematic errors in estimation of 
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Figure 11. The top panel displays the model w{9), multiplied by 9 2 from 
the top-hat bin defined by 1.0 < 2 < 1.13 with a z = 0.03(1 + z) (black 
line), <j z = 0.033(1 + z) (red dotted line) and a z = 0.027(1 + z) (blue 
dashed line). The bottom panel displays the ratio of the models with cr z = 
0.033(1 + 2) and a z = 0.027(1 + 2) to the model with a z =0.03(1 + 2). 



f(z)as(z). We note that the systematic effects are largest at scales 
greater than the BAO scale, and are therefore less likely to strongly 
affect measurements of f(z)a$(z). 



5.1 Incorrect Redshift Distributions 

One might worry that the determination of f(z)ag(z) (or any 
other cosmological parameter) will be particularly sensitive to mis- 
estimation of the true redshift distribution of galaxies within a par- 
ticular redshift bin. In order to investigate this, we alter the a z and 
also simulate the effect of a redshift dependent bias in our mod- 
elling. Previously, we assumed a z = 0.03(1 + z) for the DES 
galaxies, and we now contrast against models calculated increasing 
or decreasing this assumed error by 10%. This is about the same 
as the average uncertainty in a z stated in the DES requirements 
document, which states that the uncertainty in the dispersion in any 
photometric redshift bin of width 0.1(1 + z) should be less than 
0.003. 

The effect on the top-hat measurement of changing a z is 
shown if Fig. |TT| Its amplitude changes by ~ 10 per cent, but its 
shape remains very similar. The bottom panel of Fig.[TT] shows that 
there is a slight scale-dependency in the change in amplitude, sim- 
ilar to what C10 find (as displayed in their Fig. 6). Figs . 1 1 2| and 1 1 3 1 
display the same information as Fig. [TT| for the pair-centre mea- 
surement and the ratio of the top-hat measurement to the pair cen- 
tre measurement, respectively. The pair-centre measurement dis- 
plays almost no dependence on a z . The pair-centre measurements 
change by less than 2%, and the change in the top-hat bin is closer 
to 10%. Further, at small scales, where the redshift distortion ef- 
fect is relatively small, the ratio between the different pair-centre 
measurements differs by < 0.5 %. This implies that the pair-centre 
measurement would continue to allow for an accurate measure 
of b(z)(r&(z). However, this value of b(z)as(z) would cause the 
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Figure 12. The model pair-centre measurements for the three situations de- 
scribed in Fig. |l 1| 
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Figure 13. Same as Fig. [TT] but the ratio between the model top-hat and 
model pair-centre measurements are plotted instead, and the 8 axis is scaled 
logarithmically. 



model top-hat w(8) to differ by ~ 10 per cent from the mea- 
sured wth{6) on scales < 20/i _1 Mpc, regardless of the value of 
f(z)as(z) (since RSD do not have a significant effect on the am- 
plitude at these scales). 

The mismatch between the two measurements on small scales 
allows one to discover and fix errors in the estimation of the true 
redshift distribution, in a way that is relatively independent of RSD. 
Consider the ratio of the top-hat to pair-centre w(8): For the ratio, 
the values of b(z)a&{z) cancel each other (and one would expect 



any non-linear bias to cancel as well). This implies that even at 
small scales, where proper modelling of w(6) of galaxies requires 
knowledge of the halo occupation distribution of the galaxies, the 
ratio would remain constant. (Note, at large scales, it is the linear 
effect of RSD that causes the ratio to be non-constant.) At small- 
scales, the top-hat and pair-centre measurements should probe the 
same real-space clustering, i.e., where w(0) behaves as a power- 
law, we would expect the pair-centre and top hat measurements to 
have the same slope. Thus, we expect the ratio to be independent 
of any non-linear effects. This implies that, at small scales, any off- 
set between the measured ratio and the theoretical ratio will be due 
almost entirely to incorrect modelling of the true redshift distribu- 
tions. 

At scales where w(6) behaves as a power-law, the amplitude 
of w{9) is proportional to 



W{o z ) = / n 2 (z,a z )H(z)dz 



(28) 



Thus, one can determine how a z needs to be changed in order for 
W(a Zi a)/W(a Z: m) (where cr z ,a is the adopted mean photomet- 
ric redshift error and a ZtTn is the mean from the actual data) to 
match the offset found between the the measured and predicted 
ratios between the top-hat and pair centre measurements. When 
the errors are 10 per cent too small, the offset is 1.083, and we 
find that changing a z from 0.0027 to 0.0302 produces a match- 
ing W(a Zt a)/W(a ZtTn ). When the errors are 10 percent too large, 
the offset is 0.922, and changing a z from 0.033 to 0.0298 pro- 
duces a matching W '(a /W '(a z , m ). In both cases, a slight over- 
correction was required, but the correction would reduce the sys- 
tematic error on any f(z)a$(z) estimate to less than 1 percent. 

In practice, one will have a individual photometric redshift er- 
rors (or PDFs) for each galaxy. In this case, the correction would 
be to find the (likely constant) factor by which each error distribu- 
tion needs to widened or narrowed in order to achieve the desired 
ratio in the model. This correction is likely to be important for any 
analysis one wishes to conduct on the clustering of photometrically 
selected galaxies, and it is appealing because does not require any 
data external to the survey (unlike other proposed methods, which 
rely on smaller spectroscopic surveys, see, e.g. Newman 2008). 

We also investigate the effect of a bias on the photometric red- 
shifts (i.e., the mean photometric redshift in a bin may be higher 
or lower than the mean of the true redshift distribution). We apply 
a positive or negative photometric redshift bias of 0.01. This level 
of bias is close to 10 times higher than the 0.001(1 + z) stated in 
the DES requirements document. The effect of these biases is seen 
in Fig. 14 The model top-hat w(9) with no bias is shown by the 
black line, while red line is for zt = 0.01, and blue zt = —0.01. 
The value of the amplitudes changes only slightly — by ~1 per 
cent. This is because the main effect is a change in the median red- 
shift and thus the effective value of the linear growth rate. The ratio 
of £>(1.055)/£>(1.075) is just 1.01, thus the results we find match 
the expectation. These results suggest that the requirements of DES 
will mean that photometric redshift bias is not an issue for the pur- 
poses of measuring f(z)as{z). 



6 CONCLUSIONS 

We have described the effects of RSD on angular clustering 
measurements and outlined an effective method for measuring 
f(z)cr&(z). In order to guide these efforts, we have extended the 
techniques outlined in C10 to model the amplitude and covariance 
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Figure 14. The top panel model w(9) for a top-hat bin defined by 1.0 < 
z < 1.13. The black line assumes no photometric redshift bias, while the 
red dotted line assumes a bias 0.01, and the blue dashed line assumes a bias 
of -0.01, all for a z = 0.03(1 + z). The bottom panel displays the ratio of 
the models with bias to the model with an unbiased distribution. 



of w (9) measured using galaxies selected from a DES-like survey, 
allowing us to determine the optimal binning scheme to extract the 
most information about RSD. In particular: 

• We have shown that RSD have a measurable effect on the angular 
correlation function w(9) measured using photometric redshifts. At 
the very least, any accurate model must include their effects, or 
have shown that they are negligible. 

• We describe how using measurements w(8) in pair-centre bins 
allows b(z)ag(z) to be constrained independently of RSD, thereby 
minimizing the degeneracy between b(z)ag(z) and f(z)ag(z). 

• We have shown that f(z)ag(z) can be measured to ~ 17 x b 
per cent accuracy at z — 1 for a DES-like survey, assuming a z — 
0.03(1 + 2). 

• We confirm that our analytic predictions of the ability to constrain 
f(z)ag(z) work well as, when we measure f(z)ag(z) from mock 
catalogs generated by the MICE simulation, we obtain distributions 
that agree with our analytic predictions. 

• We show that omitting non-linear gravitational effects around 
the BAO scale in the modelling of w(6) can cause the expected 
f(z)ag(z) value one measures to be more than 15 per cent smaller 
than its true value. 

• The accuracy with which f(z)ag(z) can be measured depends 
strongly on the error on the photometric redshifts. At redshift 1, the 
expected error on f(z)ag(z) is ~ twice as large for a z — 0.05(1 + 
z) compared to a z — 0.03(1 + z). Given that the expected error 
depends strongly on the bias of the galaxies and the photometric 
redshift limits applied to the bin, we expect the optimal bin (or 
series of bins) can only be determined after the data is in hand. 

• We have shown that, after adopting a reasonable expectation 
of b(z), a series of w(9) measurements made for galaxies drawn 
from narrow (Az(l + z) < 0.02) redshift slices over the range 
0.5 < z < 1.4 should allow a DES-like survey to detect 10 per cent 
deviations from the ACDM prediction for cosmological growth of 



structure. In terms of the model f(z) = f2 m (;z) 7 , we expect aDES- 
like survey to be able to determine 7 = 0.557^Q 22- 
• We further show how the ratio can be used to correct any inac- 
curacy in the estimation of a z for a particular redshift distribution. 
This provides an attractive method for testing the photometric red- 
shifts, since it does not require any external data. 

The constraints that we predict will only improve as the red- 
shift limits and sky area increase with future surveys. Given purely 
photometric data, it is clearly still worth investigating and measur- 
ing the effects of RSD. 
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